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jrt " Abstract. We present the method of describing an evolution in quantum cosmology 

in the framework of the reduced phase space quantization of loop cosmology. We 
apply our method to the flat Friedman- Robertson- Walker model coupled to a massless 
scalar field. We identify the physical quantum Hamiltonian that is positive-definite and 
generates globally an unitary evolution of considered quantum system. We examine 
properties of expectation values of physical observables in the process of the quantum 
big bounce transition. The dispersion of evolved observables are studied for the 
Gaussian state. Calculated relative fluctuations enable an examination of the semi- 
classicality conditions and possible occurrence of the cosmic forgctfulncss. Preliminary 
*vj , estimations based on the cosmological data suggest that there was no cosmic amnesia. 

£> " Presented results are analytical, and numerical computations are only used for the 

>0 . visualization purposes. Our method may be generalized to sophisticated cosmological 
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,*~> models including the Bianchi type universes. 
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1. Introduction 

The loop quantum cosmology (LQC) method seems to be an efficient method of 
quantization of cosmological models of general relativity developed recently. Presently, 
we have two versions of this method: standard LQC (see, e.g. (TJ |2} [3] and references 
therein) and nonstandard LQC pEJ El El [7J El El \TU\ [TT]. For an extended motivation for 
developing the nonstandard LQC we recommend an appendix of our paper [6]. 

The standard LQC has been developed by several groups around the world in the 
last decade. This approach follows the Dirac program in which one first identifies the 
kinematical Hilbert space ignoring the dynamical constraints. Next, the dynamical 
constraint (choosing suitable gauges leads to a single constraint) of the theory is 
promoted into an operator acting in this space. Kernel of this operator is used to 
construct the physical Hilbert space. Similarly, observables are analyzed firstly on the 
kinematical phase space and suitable unitary transformation is constructed to map the 
kinematical results into the physical ones. 

The nonstandard LQC has been proposed recently. In the first stage of this approach 
one prepares the classical formalism for quantization: (1) Solutions to the Hamilton 
equations which satisfy the dynamical constraints are found. In the case of complicated 
system of equations, one may examine the structure of the constraint surface by the 
phase portrait methods for dynamical systems [12]. (2) Elementary Dirac observables 
on the constraint surface are determined, which define the physical phase space (PFS). 
(3) Physical observables are defined, i.e. the observables which after quantization can 
be used for predicting outcomes to be compared with observational data. Physical 
observables are introduced as functions of elementary Dirac observables and an evolution 
parameter. In the second stage, one quantizes the classical system: (1) Self-adjoint 
representations of physical observables are constructed by using the representation of the 
algebra of elementary observables. (2) Eigenvalue problems for physical observables are 
solved to get their spectra. (3) Evolution of expectation values of physical observables 
is examined. It consists in finding the new Hamiltonian of the classical theory that 
generates dynamics on the physical phase space. This new Hamiltonian is no longer 
a dynamical constraint of the classical theory. One finds a self-adjoint representation 
of this Hamiltonian. Quantum Hamiltonian is used, via Stone's theorem, to define an 
unitary operator. This operator is used to examine an evolution of the quantum system. 

Preliminary results obtained in j6j [TT] indicate that proposed method of describing 
an evolution of a quantum cosmological system is reasonable. In this paper we give 
complete presentation of the evolution. We demonstrate that it is able to reveal the 
details concerning the nature of the quantum big bounce transition. In particular, 
we analyze quantum fluctuations of physical observables in the propagation across the 
quantum bounce from the past to the future time infinities. 

In the standard LQC an evolution of cosmological system is described quite 
differently. The kernel equation for the operator constraint is used to construct, 
after some formal rearrangements, an equation interpreted as an evolution equation. 



Unfortunately, so obtained equation is usually so complicated that one can only solve 
it by combined analytical and numerical methods. This is why the preliminary 
examination of the evolution of the quantum FRW model has shown that classical 
big bang turns into quantum big bounce transition [3J, but could not say anything 
specific about the nature of the quantum bounce. In particular, an evolution of the 
dispersion effects of quantum observables was not done satisfactory. Replacing an exact 
Hamiltonian constraint of the FRW model by a simpler one have enabled making some 
analytical analysis. This simplified method for describing an evolution, called sLQC 
[15] . has shown that the cosmic amnesia for the case of semiclassical states, discovered 
earlier in analyzes of a simple cosmological toy model [H], does not occur [15] . 

Our nonstandard LQC method enables analytical studies, with an exact expression 
for dynamical constraints, of subtle quantum effects of specific cosmological model of 
the universe. In this paper, we consider the flat FRW model with a free massless scalar 
field. The choice of the model results from the fact that the dispersion effects have been 
examined so far mainly for this model. On the other hand, the model is quite simple 
and the FRW symmetry is supported by the current observational cosmology. 

In order to have our paper self-contained, we recall in Sec. II A some aspects of 
derivation of the Hamiltonian constraint of our nonstandard LQC [8]. In Sec. II B we 
introduce the notion of the physical Hamiltonian. In Sec. Ill we construct the physical 
quantum Hamiltonian and examine its spectral properties. The evolution of quantum 
FRW model is presented in Sec. IV, where we consider the dispersion of physical 
observables. We also briefly evoke the problem of time. The relative fluctuations of 
observables, for the Gaussian state, are considered as a function of time in Sec. V. We 
conclude in the last section. Appendix A includes the derivation of the formulas used 
in the section on quantum fluctuations. 

2. Classical dynamics 

2.1. Hamiltonian constraint 

The gravitational part of the classical Hamiltonian, in the Ashtekar variables (A l a , Ef), 
is the sum of the first class constraints 

H g = -^ / tfxirfCi + N a C a + NC), (1) 

where £ is the space- like part of spacetime KxS, and where Cj and C a denote the 
Gauss and the spatial diffeomorphisms constraints, respectively. For considered FRW 
model gauges are chosen in such a way that C, and C a constraints are automatically 
fulfilled. The only nontrivial part is the scalar constraint C so the Hamiltonian reads 

H ^-^iLS/vmf m ^ (2) 

where 7 is the Barbero-Immirzi parameter, and where F^ b is the curvature of SU(2) 
connection At. 



In LQC the gravitational degrees of freedom are parametrised by holonomies hi and 
fluxes Fi, which are functionals of the Ashtekar variables. These non-local functions are 
used to construct a non-perturbative theory. The holonomies and fluxes are the SU(2) 
variables satisfying the holonomy-flux algebra. In the highly symmetric spaces, like the 
FRW model considered here, the forms of these functions are simple and known. 

In particular, in the flat FRW model the flux may be parametrised by the v variable 
and the holonomy is expressed in terms of the /3 variable [6] . The variable v is a physical 
volume defined as follows 

v := / dx\dx2dx^\/ det q a b = a 3 / dxidx2dx 3 \/ det q® b =: a 3 V , (3) 

where V C £ is an elementary cell in the space with topology M 3 ; 
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fx 1 ,^ 2 ,^ 3 ) 



are Cartesian coordinates; q a b := a 2 q® b is a physical 3-metric; a is a scale 
factor; q® b dx a dx b := dx\ + dx\ + dx\ defines a fiducial 3-metric; V is a fiducial volume 
(it does not occur in final results). The /3 variable, in the limit /3 — > 0, is linked to the 
Hubble factor H = a/a via the relation (3 = ^H. 

In order to express the Hamiltonian, Eq. (jSJ), in terms of holonomies and fluxes, 
the procedure of regular izat ion has to be applied, which introduces a new scale to the 
theory, namely the parameter A. This can be understood as the length scale of the lattice 
discretization. The applied procedure of regularization and rewriting the Hamiltonian 
in terms of holonomies and fluxes is the same as known from the standard LQC (see, 
Appendix A of [8]). The obtained gravitational part of the Hamiltonian reads [3] 

ijk 

(4) 

where /in, = hihj(hi)~ l (hj)~ l is the holonomy around the square loop D^ (for more 
details see [5]), and iV is the lapse function. The elementary holonomy in the i-th 
direction reads 

h i ^cos(^y + 2smf^)r i (5) 

where r« = —\oi (o~i are the Pauli matrices). The holonomy (JSJ) is calculated in the 
fundamental representation of SU(2). The factor A is the parameter of the theory that 
may be related with the minimum area of the loop. It is expected that A ~ /pi, but its 
precise value has to be fixed observationally. 

In the model considered in this paper, the total Hamiltonian is the sum of the 
gravity Hg and free scalar field H^ := Np 2 ,/{2v) terms. The insertion of the elementary 
holonomy (J5J) into Eq. OH leads (for details, see Appendix A of [S]) to the expression^!: 

V 87r/ P1 7 2 A 2 2v) k j 

| In the rest of the paper we choose c = 1 = h and G = Ipy = 1/wpj except where otherwise noted. 



where the sign " ps " reminds that the Hamiltonian is a constraint of considered 
gravitational system. 

The Hamilton equation takes the following forml 

f = {/,# (A) }, (7) 

where 

|~ <9- <9- 9- 9- 1 9- 9- 9- 9- 
9/3 9t> 9-u 9/3 90 9p,£ 9^0 90 
For the lapse function N = 1 in Eq. (^, the time t in equation (J7J) is the coordinate 
time. Using Eqs. ([7|) and ([3]) we get the expression for the Hubble factor, H, as follows 
d 1 dv 1 sin(2A/3) . . 

"~ a ~ 3v dt ~ 7 2A ' ^ 

We have the expansion period (H > 0) for /3 E (0, ^-), and the contraction (H < 0) for 
/3 G (2A) f )■ The present classical phase of expansion corresponds to the limit /3 — > 0. 
In this limit, due to (jHJ), we get H = /3/j. This relation can be also obtained from (Q 
by shrinking the regularization parameter A to zero. 

The coordinate time t used in cosmological observations and the intrinsic time (to 
be introduced later) preferred in theoretical considerations will be related via 

at v 

Finding an explicit formula for v in terms of (see next subsection) gives more insight 
into this relationship. 

2.2. Physical Hamiltonian 

The kinematical phase space J7. of the system can be parametrized by four independent 
variables /3,v,(f), and p^. If these variables satisfy the constraint flH]), they can be used 
to parametrize the physical phase space T v which is thus three dimensional. In the 
relative dynamics one canonical variable is used to parametrize all others. Choosing 
to play such a role enables an integration of the system and finding the elementary 
Dirac observables by the method presented in [5]: 



Oi=p h (11) 

2 = - ^S^arth (cos (/3A)) , (12) 

V 12vr/pi 

/ n sin (A/3) . . 

3 = sgn(p ) v— ^, (13) 

which can be used to parametrize the phase space of the relative dynamics. Equations 
(fTTT) -(TT3l) take into account the constraint fl6]) so we have the relation [5] 

3 = kjO x , (14) 

§ / is a function on phase space. 
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where k? := AttG/3. Thus, the physical phase space of the relative dynamics is two 
dimensional. It can be parametrized by two independent elementary observables, for 
instance Oi and 0%. 

The observables are constants of motions since, by the definition of Dirac's 
observables, we have 

O k = {O kl H^} = 1 k = 1,2,3. (15) 

Thus, acting with d/dcj) on Egs. (TTTf - ffTBl leads to 

sgn(p^) A d£ 

v / 12^/piSin(A/3)^' l ° J 

= Acos(A/3)^ + sin(A/3)^. (17) 

d<p dtp 

We wish to emphasize again that (TT6l) and (TT7|) do include the constraint fl6]) . In redefined 
variables 

Q:=/3, P:=^L- V , T := -sgn(p^)v / 3^l| 1 0, (18) 

equations (TT6|) and ( Tl7|) can be rewritten formally in the form 
dP dH x 



dT " dQ 
dQ dH x 



(19) 
(20) 



dT dP ' 

where 

tf A :=-^=Psin(AQ), (21) 



plays the role of the Hamiltonian of the system devoid of the dynamical constraint!]]]. 

This way we have turned the system with constraint into the system without 
constraint. One can say that we have obtained the system in which the dynamical 
constraint has been solved. 

Let us find the solution to the system defined by ( JT9~]) - (T2~Tj) . By direct integrations, 

we find 

2 / 2 \ 

Q = — arctanexp f —j={T — T ) ) . (22) 

A WG ' 



Similarly, we get 

P = P cosh(-^(T-T )), (23) 

where Pq and To are constants of integration. 

The parameter T is the intrinsic time for the relative dynamics under considerations. 
Using fllOp . (ITS]) and ( 12 3 p we relate this time with the coordinate time t: 



t = in - 2\l^-fA sinh (-^(T - T )). (24) 



|| Proving the equivalence of Eqs. (Tl9| - ()20p with Eqs. (Tl6ll - ()17p needs using combination of (fT6|) and 

dm). 



Since Pq > (as P corresponds to the volume variable) and the r.h.s. of Eq. 
monotonically decreasing with T, one can see that the directions of t and T are quite 
the opposite. Thus, the phase of expansion in coordinate time t corresponds to the 
contraction phase in intrinsic time T and vice versa. One can give the following 
interpretation: we define the Hubble parameter, h, in terms of the intrinsic time T, 
in analogy to (j5|), to getliJ 

h ' = 3vdT = ~3 mpi cos ( A< 2)- ( 25 ) 

In contrast to the coordinate time case (Q, we have a contraction (h < 0) for 
XQ G (0,7r/2), and an expansion (h > 0) for XQ G (7r/2,7r) . Using this we can 
say again that the directions of the intrinsic time T and the coordinate time t are the 
opposite. 

It is commonly known that Hamiltonian of a physical physical system (being a 
part of some larger system) should be bounded from below, otherwise it would be 
dynamically unstable as the lowest energy state would not exist. It took much effort to 
prove that asymptotically flat spacetimes may have this property (see, e.g. the positive 
energy theorem [IE]). In the case of cosmology, the situation is quite different since the 
system, being the entire universe, is not a part of some bigger system. In what follows 
we assume that our model of the universe is an isolated system. Since we consider a 
model which is a Hamiltonian system, its total energy must be conserved. Such a system 
(with fixed spacetime geometry and topology) cannot make transition to a state with 
lower or higher energy. The only reasonable requirement is that it should have finite 
value. However, our classical Hamiltonian, defined by (21), is positive-definite since 
XQ G (0, 7r) and P > 0, which enables reasonable interpretation of the model at the 
classical level. The values of Q = and Q = tx/X can be approached by the classical 
trajectories only asymptotically. Owing to this, we postulate that the corresponding 
quantum Hamiltonian should be positive-definite too. At this stage we introduce the 
notion of the physical Hamiltonian. It is defined to be a positive-definite Hamiltonian 
that generates dynamics of the system. The property of being positive-definite is specific 
to our cosmological model and is devoid of basic importance. 

The positivity of the Hamiltonian ( T2~TT) was achieved by introducing an intrinsic 
time that has an opposite sine to the metric time. It is possible to redefine intrinsic 
time T — > —T, such that the directions of T and t will be the same. However, in such a 
case the Hamiltonian ( J2T1) would be multiplied by minus one so H\ no longer would be 
positive-definite. 

The above reasoning presents the key ideas underlying the reduced phase space 
(RPS) approach of the classical level. In what follows we present the RPS quantization. 



^f Now, wc use the Hamiltonian (|2T 



3. Quantum Hamiltonian 

In what follows we use the Hilbert space % = L 2 ([0, vr/A], dQ) so it has the scalar 
product 

(/!<?> := f V JgdQ. (26) 



The quantum Hamiltonian corresponding to the classical one fl2T|) is defined in a 
standard way to be 

mpi 



H,:~- 



X 



P sin(AQ) + sin(AQ)P 



(27) 



The classical canonical variables Q and P satisfy the algebra {Q, P} = 1. Choosing 
the Schrodinger representations for these variables 

. d 



Q<P(Q) := Q(f>(Q), P<j>(Q) := 



'dQ 



(28) 



where <fi G "H, gives formally [Q, P] = il. 

An explicit form of the operator (T2~T|) is easily found to be 



H x ^ = -j±= (2sm(\Q)^L + Acos(AQ) ) V, 



(29) 



where ip G D C "H, and where .D is some dense subspace of "H. In what follows we wish 
to determine D which is the domain of self-adjointness of the operator H\. 



3.1. Eigenproblem for the Hamiltonian 

The eigenequation for the operator H\ reads 

i / d^ \ 

- 3 - 7 =(2sm(AQ)- + Acos(AQ)*)= M . 

The solution of the above equation is given by 

where \l/o is a constant of integration. Let us calculate 



(Ve\Ve>) = |*c 



-A exphiy/G(E' -£)ln|tan(^)|} 



dQ 



sin(AQ) 



Defining 



G 



x :- 



In 



tan 



\Q 



once can rewrite the above integral into the form 

2 r+ao 



[*e\*e>) = l*o| 2 -^= / exp {i(E' - E)x} dx 
XyG j-oo 



l^r 



XV G 



2tt6(E'-E). 



(30) 



(31) 



(32) 



(33) 



(34) 



By choosing 



\/^P ( 35 ) 



we get 

(36) 



MQ) = iQ\**) = V^Jgy-p [liVGEln |tan ( 



AQ 
2 

The states ^ E given by equation ( l36|) satisfies the orthonormality condition in the form 
(*s|^) = 8{E' - E). (37) 

Based on this, for the eigenstates, we obtain 

{^ E \H^ E ,) - (H>$ E \$ EI ) = (E' - E)5(E' -E) = 0, (38) 

which means that the Hamiltonian is symmetric on the space of the eigenstates. 

3.2. Symmetricity 

We define the domain of H x as follows 

D(H X ) := span{<^, k G Z}, (39) 

where 

/oo 
c k (E)* E (Q) dE, c k G C °°(R). (40) 

-oo 

It is clear that D(H X ) C £ 2 ([0, vr/A], dQ) = % is a dense subspace of "H, and an action 
of the unbounded operator H x does not lead outside of D(H X ). 
For any tp k and (pi we have 



{tp k \H\tpi} - {H x <p h \<pi) 

-+oo P+OO 



/+oo /"+oo 
/ c k (E) Cl (E') (E' - E)5(E - E') dEdE' = 0. 
-OO J —OO V 

=0 

Therefore, the operator .£/>, is symmetric on D(H X ). 



(41) 



5.5. Self-adjointness 

To examine the self-adjointness of the unbounded operator i^Aj we first identify the 
deficiency subspaces, /C±, of this operator [18] 

£± := {g± G D(^) I (g±\(H x ±U)<p,) = 0, \f V e D(H X }, (42) 

where D(^) := {/ G L 2 ([0, tt/A], dQ) | 3!/* (/*|<7> = (f\H x g), \/g G £>(#a)}, and 



where 



/oo 
b(E)* E (Q)dE, 6gC °°(R). (43) 

-oo 

For each y9 fc G D(H X ), defined by (1401) . we have 
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0={g ± \{H x ±iK)<p k ) 

P7T/X pco pco 

= dQ dE l dE 2 b* ± (E 2 )^* E2 {Q){H x ±H)c k (E 1 )^ El {Q) 

JO J — oo J — oo 

/oo pco pir/X 

dE x / dE 2 b* ± (E 2 )c k (E 1 )(E 1 ±i) V* E2 (Q)V El (Q)dQ 

oo J —oo Jo 

/oo 
dE(E±i)b* ± (E)c k (E) =► 6 + = = &_ . (44) 

-oo 

Thus, the deficiency indices n± := dim[K,±] of H\ satisfy the relation: n + = = n_, 
which proves that the operator H x is essentially self-adjoint on D(H\). One can argue 
that the whole spectrum of H\ belongs to the real line. 

3.4- Physical Hamiltonian 

The classical physical Hamiltonian H\ is positive-definite. The corresponding self- 
adjoint operator H\ has however eigenvalues E G K. We therefore introduce the 
quantum physical Hamiltonian H by requiring that it has only positive eigenvalues. 
It is defined as follows 

W\V E ) ■= \E\\Ve), E>0, (45) 

where \^e) is the eigenvector of the Hamiltonian H corresponding to the eigenvalue 
E. It is clear that the spectrum of the operator HI is doubly degenerate. For any state 
$ = J+™ c(E)^ E dE G D(H) = D(H X ), where c G Cg°(R), we get 



-oo 

™+oo 



/-t-oo 
dE\c{E)\ 2 \E\ >0. (46) 

■oo 



4. Evolution of quantum system 



Making use of the Stone theorem [18] , we define the unitary operator of an evolution as 
follows 

U(s) := e^, (47) 

where s G M is a 'time' parameter. The state at any moment of time can be found as 
follows 

|*( 5 )) = U(s)\V(0)) = e- isfi |*(0)). (48) 

The minus sign in (J4T1) . of the exponential function of an operator, is essential as only 
in this case an infinitesimal version of U(s) acting on |\l/(s)) leads to the Schrodinger 
equation. 

Let us consider a superposition of the Hamiltonian eigenstates 

|#(0)> = / dEc(E)\V E ) (49) 
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at s = 0. Then evolution of this state is given by 

/+oo 
dEc(E)e~ isM \^ E ) 
■oo 
/+oo 
dEc(E)e- ia \ E \\V E ). (50) 

■oo 

Normalization of this state is given by the following condition 

/+oo f+cx> 
/ dEdE'4E)c(E')e- is ^-\ E 'V (V E \V E ,) 

=5{E'-E) 

/+oo 
dE\c{E)\ 2 = l. (51) 

-oo 

Now, let us assume that the superposition (1301) has a form of the Gaussian packet with 
the profile defined to be 



2a 



1/4 



c(E):=l^\ e -(E-E )^ (52) 

centered around E with the dispersion parametrized by a. The normalization factor 
(2a /tt) ' is due to the condition (l5"Tj) . In the rest of this paper we study properties of 
this state only. Assuming that Eo^/a ^> 1, one can approximate 

*(Q, s) = (Q|*(s)) ( — J / dEe-^-^-^^tQ) 

2aV / A ^ /" +0 ° j R - a (E-E r-isE iEx ,,o^ 

J V 4 ^ MAC) 7-00 ' ( } 

where a; := |\/Gln |tan (-^) |. The approximation is valid because, for E oy fa 3> 1, 
contribution from the negative energies to integral ( 1531) is negligible. Therefore, |i£| can 
be replaced by £7. Calculating the integral (!53|) we get 



M / - ? i=- X e-^e^^-), (54) 

V 87ra y sin (AQ) 

where a := a/G. The probability density is easily found to be 



|*(Q,s)' : 2 * 



IH tan (f)l -&' 



(55) 
A v / 8™sin(Ag) 

It is worth to emphasis that the state (154"]) is normalizable so the probabilistic 
interpretation can be applied. Thus, the function \^?(Q, s)\ 2 gives the probability density 
of finding universe with given Q at the particular moment of time s. In the case 
A — > 0, the state becomes non-normalizable and the probabilistic interpretation cannot 
be applied. Finding a normalizable state is a serious problem in most of quantum 
cosmologies different from LQC (see, e.g. [TT]). 
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4-1. Evolution of observable Q 

One can now investigate the mean value of the Q operator in the |^(s)) state. It is 
clear that Q is a symmetric bounded operator on |\I/(s)). We find 

(Q) := (*(s)|Q|*(s)> 



arctan < exp 



2 [y 



h'l 



One can show that 

lim X(Q}(s) = 7r, lim \(Q}(s) = 

s— >+oo s— »— oo 

which agrees with the classical limits. We also find 

(*( s )|g 2 |v|,( s )) = 



exp 



h? > dy - 



V2 



ira 



A 



+ 00 



arctan < exp 



2 [y 



tpi 



■^v 2 



dy. 



Thus, the dispersion of Q in the state |^(s)) reado 



AQ := V(*00IW(*)> - ms)\QMsW. 



(56) 



(57) 



(5* 



(59) 



In Fig. [T] we show (Q) (thick blue line) and compare it with the classical solution 
(dashed red line) 



2 / 2 

Q =- — arctan exp f —= (T — To 



with s = T — Tq. The shadowed region represents the dispersion of our state, and it is 




Figure 1. The thick blue line represents the mean value (Q) for a = 0.1. The dashed 
red line is the classical solution Q = j arctan exp {2s/ipi}. The shadowed region is 
constrained by (Q) ± AQ. 



The integrals (|56|) and (|58jl have been determined numerically. 
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constrained by (Q) ± AQ. One can notice that, close to the bounce at s = (where 
XQ = 7r/2), the dispersion of the state is significant, while it fast decreases ass-> ±00. 
Therefore, far from the bounce, the state converge to the classical solution_|. 



4-2. Evolution of observable P 

The volume operator, due to ({TBI , reads 

v = 47r/p 1 7P. 
The P operator is unbounded, but symmetric on the state |^/(s)) defined by Eq. 



(61) 
51): 



(*( S )|Ptf( S )) 



A 



lim cosh 



\Z8na a->+oc 

+<P¥( S )|*(s)> = <P*(s)|tf(s)>, 

1 



G 



.r 



2c — e 



2a 



where we used the relation cosh ( -j=x 



sin(AQ) ' 

It is not difficult to derive the following: 

(P) := (¥( S )|P|¥(*)> 

AP f + °° ,( 2 \ 
= = / cosh ( —j=x ] exp 

V87ra 

1 



G 



la 



[x 



dx 



AP /pie 25 cosh 



2s 



Therefore, 



(v) = 4t4 i7 (P) = 27r4 l7 AP e 25 cosh 

The minimum allowed volume in the \& state is 
(v) min = 2nl 3 Pll \E e 2 \ 

The corresponding classical solution reads 

/ 2 
v — Vn cosh ( —= (T — T 



2.s 



G 

Thus, the functional forms of (v) and v coincide. We also find 



;*( s )|p 2 |*( s) ) 



A 2 

32 



4 + 4P; 



cv 



+ cosh 



4.s 



,85 



GJ \a 



- + 4 + 4P 2 



where Po := Po/Ppi- The dispersion of P in the state |^(s)) is found to be 

AP:= V / (^(s)|P 2 |^(s)) - ((^(s)|P|^(s))) 2 



A 



2(4 



4 + 4E7jJ + 2 cosh 



is 



,8a 



- + 4 + 4P 2 



(62) 



(63) 

(64) 
(65) 
(66) 



(67) 



16P 2 e 4fi cosh 2 



2.s 



1/2 



* One can notice some little discrepancy between the classical solution and the mean value (Q) 
is however due to the particular form of the state ([52]) that has been chosen. 



(68) 
This 
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We visualize the dispersion of P as a function of time in Fig. [2J 



lor 
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Figure 2. The thick blue line represents the mean value (P) for a 
The shadowed region is constrained by (P) ± AP. 



0.1 and E = 10. 



We can see that both the volume and its dispersion grow quite fast away (exponentially) 
from the big bounce regiodj. This can be understood in the context of Heisenberg's 
relation. Namely, while dispersion AQ tends to zero, for s — > ±oo, the dispersion AP 
grows appropriately to fulfill uncertainty relation AQAP > 1/2. We shall study this 
issue in more details in the next subsection. 



4-3. The Heisenberg uncertainty relation 

The standard probabilistic interpretation of quantum mechanics cannot be applied to a 
cosmological system for a number of reasons. For instance, there is only one Universe 
and there is no observer outside the Universe. Thus, the process of measurement in 
quantum cosmology may differ from that known from the Copenhagen interpretation of 
quantum mechanics. Instead of complaining at the interpretation problems, it makes 
sense verification if some fundamental relations underlying quantum mechanics are 
satisfied. The best known is the Heisenberg uncertainty relation. What is its status 
in our cosmological setup? Is it satisfied during the quantum big bounce transition? 

Our two canonical variables Q and P do not commute: [Q, P] = il. Since both 
operators are symmetric on the space of states |\P(s)), they should satisfy algebraically 
the Heisenberg relation: 

AQ AP > 1/2. (69) 

In Fig. |3]we show the evolution of AQAP for the model considered in this paper. 
We can see that the Heisenberg relation ( )69|) is perfectly satisfied during the entire 

(t The P-axis uses the logarithmic scale. 
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Figure 3. Evolution of the product of dispersions AQAP. Heisenberg's uncertainty 
relation AQAP > 1/2 is preserved during the evolution. 



evolution, for any values of the parameters a and E of the state. However, one cannot 
ascribe probabilistic interpretation to this relation. It is so because the representation 
[Q, P] = zl cannot be self adjoint owing to the fact (shown earlier) that Q and P are 
bounded and unbounded operators, respectively [19] . 

The function AQAP is symmetric with respect to the bounce and reaches its 
minimal value at the bounce. This is rather surprising result. Namely, one would 
expect that the transition point (s = 0) is the most quantum part of the evolution, 
while in the limits s — > ±oo one should get the most classical evolution. The situation 
is however quite the opposite. The transition point (big bounce) is the least quantum 
part of the evolution! In turn, in the limits s — > ±oo the product of dispersions AQAP 
is saturated: 

lim AQAP = 

s— >±oo 



^ oV ^e«-l + e«(-|_ + -L)4 (70) 

which shows that the Gaussian packet we consider is always quantum. We can also 
see that the relation (I7CT|) does not depend on the parameter A, which is a remarkable 
feature of our quantization schem affl 

Let us investigate in more details the value of AQAP at the bounce, AQAP^. In 
Fig. H] we show AQAP\^ as a function of a for E Q = 1, 10 and 100. As we can see, 
the boundary value AQAP\^ = 1/2 is never crossed and it is approached for a — > 0. 
Therefore, the more sharply peaked the state is, the smaller value of AQAP is at the 
bounce. For any value of a, the smaller value of E is the boundary AQAP\b = 1/2 

ff The parameter A appears in the formalism as the result of approximating the curvature of connection 
by holonomics around small loops with length A. It is a free parameter of the nonstandard LQC and 
fixed parameter of the standard LQC, so it is of basic importance. 
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Figure 4. The value of AQAP at the bounce (s = 0) as a function of the parameters 
of state. 



is easier approached. In order to understand this dependence better, let us investigate 
separately AQ|b and AP|b- We show these dispersions, as a function of a, in the left 
and right panels of Fig. [51 respectively. The function AQ|b grows monotonically with 



0.5- 




a) 



io ■- 



b) 



Figure 5. a) Dispersion XAQ at the bounce (s = 0) as a function of the parameter 
a. b) Dispersion AP/X at the bounce (s = 0) as a function of the parameter a for 
E = 1, 10 and 100. 



increase of a and is independent on E . The dependence on a is more complex for AP|b- 
Namely for any Eq there always exists some a at which the function AP|b(<5) takes the 
minimum. The smallest possible value of AP|b ~ 1.0686A is reached for Eq — > and 
a « 0.11797. 

In summary, our quantum cosmological setup is devoid of complete standard 
probabilistic interpretation, but satisfies formally the most basic relationships of the 
standard quantum mechanics. 
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4-4- Problem of time 

An evolution of classical variables Q and P, as presented by Eqs. ( |22l) and ( l23l) . is 
parametrized by a free massless scalar field 0, due to ( fl8l) . that is a monotonic function so 
it can play the role of an internal clock [5] . The expectation values of the corresponding 
quantum operators Q and P, defined by ( 1471) and ( 1631) . are parametrized by s £ IR, owing 
to (j4"Tj) . Are the evolution variables and s quite independent? An important difference 
between them is that they label an evolution of the system at different levels: classical 
and quantum, respectively. We postulate that these two variables are related linearly: 
s = ai<p + a2, where ai,a 2 £ M. It means that neither <fi nor s belong to the physical 
phase space. This seems to be the specific feature of our reduced phase space (RPS) 
method. It has been already proposed in [7j treating the <p variable as an evolution 
parameter of both classical and quantum dynamics. Such an interpretation is supported 
by the plots of Fig. [T]and Fig. [21 where the same abscissa is used to label both classical 
and quantum evolution of presented functions. The plots of classical and corresponding 
quantum functions practically coincide. Such an agreement suggests that our postulate 
concerning time variable is reasonable. 

5. Relative fluctuations 

In this section we study the relative fluctuations of the quantum observables O in 
the state |\l/). We consider three observables: HI, Q and P. The relative fluctuation 
AO/(0) is a measure of the semi-classicality of a quantum state. We say that |\&} 
is semiclassical if AO/(0) <C 1, and quantum if AO/(0) ~ 1. It is clear that the 
simiclassicality notion is not at all defined uniquely. We apply such a definition of 
semiclassicality because we wish to be able to make comparison of our analyzes with 
the available published results [HI [151 1201 I2E 122] • In the future we shall try to apply 
a more sophisticated definition: If the uncertainty in an observable is less than the 
observational precision, the state is semiclassical with respect to that observable. 

In what follows, we also consider the function Dq characterizing the asymptotic 
aspects of relative fluctuations with respect to the bounce, defined to be [15] 



Dq := lim 

s— >oo 



AO(-s 



(o)(- 





(71) 



5.1. Relative fluctuations ofW 



Let us start from deriving expectation value of HI in the state |\l/(s)), we find 

(HI) := (tt( s )|H|tt(s)> = E erf (v^Eq) + ^=^, (72) 

where erf(x) := -7= J x e _t dt is the error function. One can see that for \faE ~^> 1, the 
above expression simplifies to (HI) ~ E . In order to find the dispersion of EI, we also 
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need: 



[y(s)\M 2 \y(s)) = ± + El 



(73) 



Based on the above, we determine 



AH := V(*(s)|H 2 |^(s)) - ((^(s)|H|#(s))) 2 
1 



Act 



2E e 



Ei 



1 - erf 2 ( v^Eq 



-IolE}, 



V2 



ira 



erf V2aE 



2-Ka 



1/2 



(74) 
(75) 



We can see that for y/aEo ^> 1, dispersion of EI simplifies to AH ~ 1/yAa. Therefore, 
for yfaE ^> 1, the relative dispersion 

AH 1 



2^„ « »■ < 76 ' 

It is worth to note that condition ^ s /aEo ^> 1 was also used while performing integration 
( 153]) . One can see now that approximation based on this condition is justified by 
the restriction imposed on the relative fluctuations of the Hamiltonian H. The 
semiclassicality requires \/q.Eq 3> 1. The relative fluctuations of H are constant in 
time and therefore symmetric with respect to the bounce so finally we get: D-a = 0. 



5.2. Relative fluctuations of P 

Relative fluctuations of P can be expressed as follows 



AP 



\ 



4 - 4 + 4£ 2 ) + cosh (^)e te (lH AE? 



8E$e A& cosh 2 (^ 



1. 



(77) 



In the left part of Fig. |6] we plot this relation for the different values of a and E$. 
The relative fluctuations of P are symmetric with respect to the bounce and reach the 



a) 



S = 0.15, £„ = 10 



0.8 5= 0.1, E„ = 10 





-1 -1 U 1 1 i 

»M b) 5 

Figure 6. a) Relative fluctuations ^£ for different values of a and -Eo- b) Relative 
fluctuations j£ at the bounce (s — 0) as a function of a for i?o = 1, 10 and 100. 
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minimal value at the transition point between the contracting and expanding phases. 
The symmetry s —)■ —s directly implies that 



Dt 



lim 

s— >oo 



AP(-s) 



AP(s) 



0. 



(P)(s)J \(P)(s 
Thus, at infinite past and future times, the relative fluctuations of P are the same: 



(78 



AP 



(P) 



lim —^~ 

s->±oo /p\ 




1 + e 



io 



4Ega 



El 



(79) 



The function j£- 



(a) has the minimum for any Eq, which is located at 

-1 + J5 + 4E% 



a, 



8(1 + ED 



Therefore, having E , one can always minimize j£ 



a, 



0) 
If 



(a) by choosing a 

Eq ~^> 1, the expression (1501) can be approximated by <S min ~ 1/(4P ). 

Let us now consider the fluctuations of P at the bounce, which can be expressed 
as follows 



AP 



(P) 



\ 



8P 2 e 4fi 



1. 



We plot this function in the right part of Fig. [0J for fixed values of E = 1, 10 and 
100. As we can see, for any value of E , there is some a at which fluctuations take the 
minimum. These minimal value fluctuations decrease with the increase of Eq. In the 
limit Eq — >■ oo, the relative fluctuations at the bounce are given by 

AP 



lim 

So^oo (P) 



V2sinh(2a). 



(82) 



Therefore, the fluctuations P at the bounce go to zero for a — > and E — > oo. 



5. 3. Relative fluctuations of Q 

We have found (see, Appendix A) that after the bounce the relative fluctuations are 
decreasing and go to zero, so we have 
AQ(s 



lim 

-h-°° (Q){s) 



0. 



3) 



Therefore, the state becomes a semiclassical one. However, it may not be the case in the 
far past for large enough a. While moving backward in time the relative fluctuations 
saturate. This saturated value can be found (see, Appendix B) by calculating the 
integrals: {^(s)\Q\^f(s)) and (^(s)\Q 2 \^(s)) . One gets 



AQ 

(Q) 



lim 



AQ(s) 
(Q)(s) 



\/e 4 ° - 1, 



n) 
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where a > 0. Thus, before the bounce the state may be a quantum one if the value a 
is sufficiently large. 

In Fig. 0we present the plot of AQ/(Q) for different values of a. The relative 
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Figure 7. Relative fluctuations — ^ for different values of a 



fluctuations are not symmetric across the bounce (s = 0). For s — > +oo the relative 
fluctuations converge to zero. For s — > — oo they saturate giving y/e 4& — 1. The 
difference in the asymptotic values of the relative fluctuations is found to be 

D Q = e 45 - 1. (85) 

In order o interpret these results it is crucial to recall that directions of the parameter of 
time s and coordinate time t are opposite. Therefore, the positive values of s correspond 
to contraction while negative to expansion. Therefore, the relative fluctuations of Q grow 
from the contracting to the expanding phase. So, the universe becomes more quantum 
with the increase of time t. The fluctuations tends to zero for t — > — oo. Therefore the 
universe started its evolution from a sharply peaked state. Let us call it a semiclassical 
state. Only if the value of a is sufficiently small the universe will keep its semiclassicality 
during the whole evolution. 



5.4- Semiclassicality 

One can find that the maximal relative fluctuations of P and Q are related as follows 



AP 



(P) 



AQ 



<Q> 



-, i" 



i 



AEI& El 



Thus, we have 



AP 

W) 



> 



AQ 

(Q) 



(86) 
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The equality is obtained in the limit E — > oo, which leads to the equation: 



AP 



(P) 



Ve 4a - 1 



AQ 



(Q) 



It is clear that the semiclassicality condition 
AP 



(P) 



<1, 



implicates, due to (1ST)) , that we have 

aq 



(Q) 



<i. 



(89) 



(90) 



Therefore, the semiclassicality imposed on P guaranties the semiclassicality for Q. We 
can see that there is no cosmic forgetfulness if the condition (IS9)) is satisfied. 

The relation between maximal fluctuations j£ and the fluctuations at the 



bounce -£- 



is given by 

/ \ 2 

f AP 



\(P) 



tanh(4a) + -= 



,1q 



+ 



AP 



P 2 cosh(4a) cosh(4a) \(p) 



which leads to 



AP 



Jla 



> 



AP 



' v / cosh(45) (P) 
Therefore, if the maximal fluctuations of P are constrained: 
AP 



then we have 



(P) 



AP 



<1, 



(P) 



< e 26 v / cosh(4a) < 1, 



(91) 



(92) 



(93) 



(94) 



which proves (owing to -S <C 1, which implies a <C 1) that the bounce is 

W) max 

semiclassical as well. 

<C 1. Can this 



We have shown that condition -^ 



<C 1 implies ^ 



implication be true also in the opposite direction? The condition jM 



(Q) 



< 1 it 



equivalent, due to ( 1S4)) . to the restriction a < 1. By taking this into account, the 
maximal relative fluctuations of P can be expressed as follows 



AP 



(P) 



2^ 
Po 



0(a 



^/ 2 ) 



(95) 



2 y/aE 

The first term in this expansion grows rapidly with decrease of a. Therefore, unless 
the value of E is not sufficiently large, the condition j£ ^C 1 is not fulfilled. This 

\ ' max 



22 



condition is fulfilled if En ^> l/\/~a, but this is exactly requirement of the semiclassicality 
of the relative fluctuations of HI. Therefore, if the condition ( ITS"]) is fulfilled, we have the 
equivalence: 



AP 



(P) 



<1 



AQ 



(Q) 



<1 



(96) 



To complete our considerations, we show that the saturated value of AQAP can 
be easily expressed in terms of maximal relative fluctuations of Q and P: 



lim AQAP = -e 4& E ^S- 

«->±°° 2 (Q) 



AP 



(P) 



1 
> -. 

~ 2 



(97) 



Owing to Heisenberg's uncertainty relation, the above expression leads to the following 
constraint 



AQ 



(Q) 



AP 



(P) 



-4a 



> 



En 



This constraint, together with inequality ( 18 7p . results gives 



AP 



(P) 



-2a 



> 



E 



(98) 



(99) 



Thus, the samiclassicality condition ( 1891) leads to the following constraint: 

E > e~ 2& . (100) 



5.5. Forgetfulness 

Let us try to answer the question: Was the Universe quantum or semiclassical before the 
big bounce? This question is related to the problem of cosmic forgetfulness discussed 
recently in papers [T^] [T5] 120] |2T] 122] . If the Universe was quantum before the big bounce 
and the bounce transition turned it into semiclassical one, we can talk about a sort of 
cosmic amnesia. In this case the complete information about the Universe before the 
bounce cannot be obtained from the observational data after the bounce. 

The above question has been addressed so far mainly by an examination of the 
relative fluctuations of the volume observable v (proportional to our P observable) 
[2"T]|2"2"]. However, as we have shown, this type of fluctuation is symmetric with respect 
to the bounce. Therefore, the constraint on relative fluctuation at some time +\s\ 
results with the same constraint on the fluctuations at the time —\s\. As one can see on 
Fig. El the fluctuations of P saturate quickly outside the neighborhood of the bounce 
(within a few Planck's times). Thus, present cosmic observation of the semiclassicality, 
tst = fly *C 1, of the Universe would guarantee its semiclassicality in the distant past 
before the bounce. However, the situation is that this type of relative fluctuation cannot 
be 'measured'. The reason is that one does not know, first of all, how to measure the (v) 
quantity. It could be possible to measure the volume, to some extent, if the Universe was 
curved and the curvature term was measured in astronomical observations. At present, 
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there is however no indication for such a contribution. Therefore, the fluctuations of v 
are not measurable so the cosmic forgetfulness cannot be examined by using the relative 
fluctuations 4r. 

(v) 

What about the Q observable? The variable Q is directly related to the 
expansion rate, i.e. the Hubble factor fl9]), which is a quantity that can be determined 
observationally. Thus, the value of (Q) can be measured. Also the observational 
uncertainty of the Hubble factor can be used to constrain AQ. Therefore, relative 
fluctuations of Q can be constrained observationally! The present value of the Hubble 
factor is H = 70.2 ± 1.4 km s" 1 Mpc" 1 [23], therefore ^§^- « 0.02. As we have shown 
in Sec. II, in the classical limit Q = / ~fH. Thus, we propose to consider the constraint: 
-S- < 0.02. This is because the relative quantum fluctuations cannot be greater than 
the relative uncertainty of measurement. 



Due to the relation (124")) . the directions of the intrinsic time T and the coordinate 
time t are the opposite. Therefore, relative fluctuations of Q at T — > — oo correspond to 
the limit t — > +oo. Thus, the relative fluctuations of Q grow in the coordinate time t and 
saturate at y/e i& — 1. The model we consider is applicable to the evolution in vicinity of 
the Planck epoch. However, if we assume that the relative fluctuations behave similarly 
threafter, the present restriction — Q < 0.02 can be used to constraint the model. The 

condition \/e 4fi — 1 < 0.02 translates into the condition a < 10~ 4 . If the Gaussian state 
with such value of its parameter can be treated as a semiclassical one, we can say that 
the amnesia does not occur. In such a case the contraction and the bounce phases are 
semiclassical, so we have the cosmic recall. 

The above constrains are quite preliminary. There exist the possibility to put 
more robust constraints based on the phase of inflation and observations of the Cosmic 
Microwave Background Radiation. However, for this purpose the model has to be 
generalized by taking into account potential of the scalar field. 

6. Summary and Conclusions 

In the Hamiltonian formulation of general relativity, GR, the total Hamiltonian is a 
linear combination of the constraints (see, e.g. Eq. (Q])) so it cannot play the role of 
the generator of an evolution of a gravitational system. On the other hand, the GR 
system evolves according to the Einstein equations. Can one overcome this difficulty of 
the Hamiltonian formulation? There are two ways of dealing with this problem: 

(i) One eliminates the time variable in favor of any canonical variable by some formal 
trick carried out on Hamilton's equations (see, e.g. [5] for more details), which 
leads to the so called relative dynamics commonly used in LQC. In this procedure 
the constraints are used to define the physical phase space. However, the evolution 
is poorly defined since one gets the dependance of canonical variables in terms of 
any specific variable so the dependance on time may become deeply hidden. 

(ii) Using the constraints, one expresses the specific canonical variable (of Hamilton's 
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equations) in terms of other variables. Elementary Dirac observables are constants 
of motion and include the dynamical constraint. The new Hamilton's equations 
(including the constraint) are used for finding the Hamiltonian that generates 
dynamics on the physical phase space. This new Hamiltonian is no longer a 
dynamical constraint of the classical theory, but the generator of dynamics, formally 
free of dynamical constraints. This approach restores the notion of an evolution of 
a classical system. 

One should remember that the above considerations on the evolution parameter 
apply first of all to the situation with 'matter' field, like the scalar field, that can be 
used to define this parameter. It would be interesting to extend these ideas to the case 
with no sources. The simplest nontrivial example is the Kasner model. 

When we wish to quantize a Hamiltonian system with constraints, we may apply 
the Dirac or the RPS quantization methods. Both methods are plagued by numerous 
ambiguities. Our analyzes are devoid of the need of any restriction to some superselection 
sectors that naturally arise in the standard LQC quantization. The latter need not be 
a drawback of the method as it may leave some interesting imprints in the physical 
predictions [21]. The former one, seems to be less complicated than the Dirac method, 
and offers an analytical insight into physical aspect of considered model. Presented 
results demonstrate, to some extent, that our quantization scheme leads to a quantum 
cosmological system with general properties of a quantum systems we are dealing with 
in terrestrial laboratories. 

It seems that one can apply our method to much more complex cosmological models 
than the FRW universe. Recently, we have managed to quantize the Bianchi I model 
[T0| [TTj . The case of the Biachi II model can be treated by analogy. In summary, we 
suggest that quantization of cosmological systems in terms of the RPS method is much 
more efficient and unique than Dirac's method. However, since quantum cosmology 
'experimental' data are not available yet, the best strategy seems to be applying both 
methods to compare the results. An agreement of the results would prove that the 
procedure of quantization was correct. 

The loop quantization methods applied to simple cosmological models teach us that 
approximating the curvature of connection by holonomies around small loops enable 
replacing classical sinularities by quantum bounces [IJ|2l[3j|4j[5j|6l[7llEl[9l [101 El- 
To get some information about the nature of a bounce, one examines propagation 
quantum states across the bounce [151 EH 1211 122]- Such method is similar, to 
some extent, to the method used, for instance, in nuclear physics where one scatters 
a particle against an atomic nucleus to get information on the structure of the latter. 
In papers [HI |2Q] one considers solvable toy models (motivated by LQC) to argue that 
a quantum state before the bounce may become decohered at the bounce and become 
semiclassical afterwards. Applying the sLQC prescription authors of [151 121] claim that 
there is no cosmological amnesia at all: suitable semiclassical state before the bounce 
keeps being semiclassical after the bounce. Authors of [22] give strong support to this 
result by applying various analytical and numerical methods, within the standard LQC, 
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to general forms of semiclassical states. They identify the condition under which one has 
the preservation of the semiclassicality across the bounce. However, they have mainly 
examined the volume observable which is little useful for testing the cosmic amnesia. 

In our paper we have considered the transition of the Gaussian wave packet across 
the bounce of the quantum FRW universe within an exact framework. It results 
from our studies that the Q observable is the proper quantity to study the cosmic 
amnesia. It is because relative fluctuations of Q can be observationally constrained. 
Moreover, the semiclassicality condition imposed in the expanding phase restricts also 
the quantum fluctuations in the contracting phase. The preliminary observational 
constraint —Q < 0.02 indicates that the semiclassicality condition, as defined earlier, 
may be fulfilled. Owing to this, one can infer that there was a cosmic amnesia or there 
was a recall depending on what we mean by a semiclassical state. Our results support 
the prediction of the standard LQC [15], EH, [22] . We suggest to repeat the calculations 
with the variety of states different from the Gaussian type states to verify our results. 

On the other hand, the LQC results obtained for the FRW type models cannot be 
probably used successfully to describe the Universe. The very high symmetry of space 
specific to the FRW model is probably unrealistic near the cosmological singularity. We 
suggest that the real nature of the bounce may become known only after we quantize the 
Belinskii-Khalatnikov-Lifshitz (BKL) scenario [25 ] 126 ] 127], which concerns the generic 
cosmological singularity. Quantization of simple cosmological models carried out during 
the last decade may be treated as warming up before meeting this challenge. 
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Appendix A. Quantum asymptotics of Q. 

It this appendix we study dispersions and relative fluctuations of observable O in the 
limits s — ¥ ±oo. We show these limits can be found analytically by performing suitable 
expansions of integrals in expressions ( 156]) and (158]) . Based on this, we derive equations 
(183]) and (El]). 

Appendix A.l. The case s — > +oo 

—2-2- 

Let us introduce e := e *pi, which tends to zero for s — > +oo. Based on this, one can 
perform Taylor expansion with respect to e, as follows 



e 2y I vr 



arctan<^ — } = - - e~ 2y e + 0(e s ). (A.l) 
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This expansion applied to equation (1551) gives 

A A 
By squaring expansion (lA.ip . we obtain 

arctan 2 j-^j = (iy -ne-^e + e-^e' + Oie 3 ). (A.3) 

This expansion, applied in equation (1581) . leads to 

n (l) e 2& e+(^] e 4 V + (9(e 3 ). (A.4) 



,A/ V^/ VA 

Based on ([Oj) and f[Q|) . we find 



(Q 2 ) ~ ((Q)f = 0) 2 e 45 (e 4fi - l) e 2 + 0{e% (A.5) 



which leads to the following expression for the dispersion of Q: 



A<? = \J{Q 2 ) ~ ({Q)) 2 = ^e 2 Ve 4i -l + 0(e 3 ). (A.6) 



Using (IA.2P and (IA.6[) we get 

~ 2 



lim — *— = lim 

e ^° (Q) e ^° 



e-e 2 Ve 43 - 1 + 0(e 3 
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(A.7) 



which proofs equation (1831) . 

Appendix A. 2. The case s — >■ — oo 

Let us introduce the variable e := e 'pi, which tends to zero for s — > — oo. It is worth 
to stress that the parameter e introduced here differs from that used in the previous 
section. We perform the Taylor expansion 

arctan {e 2y e} = e 2y e + C(e 3 ). (A.8) 

This expansion, applied in equation (IB'81) . leads to 

1 2 [ +oc f 1 • 1 

(Q) = /^-t T / arctan {e 2y e} exp <^ -— y 2 \ dy 
\j2ixa A J_ 00 { 2a ) 

= ^e 2 "e + C(e 3 ). (A.9) 

A 

By squaring ( 1A.8I) . and applying it to equation ( )58l) . we find 

{Q2) - vk 2 (!) 2 C MCtm2 { A} exp {-^1 d!/ 

I) e 8( V + 0(e 4 ). (A.10) 



Expansions (1A.9|) and ( lA.lOj) lead to 

(Q 2 )-((Q)) 2 =(f) 2 e 4fi (e 4fi -l)e 2 + 0(6 4 ). 
Based on this, dispersion of Q in the limit s — > — oo is given by 



\l (Q 2 ) - ((Q)) 2 

With use of ljO)l and fjXl2|l we find that 
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ef e^Ve 4 " - 



, A <0 , e§- 2( ^/^ 
lim — — = lim 

«~0 (Q) e ->0 



(9(e 2 



2 2a 



C(e 3 



lim 

e-S>0 



Ve 4 * - 1 + C(e 2 



V^ 



la 
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(A.ll) 

(A.12) 

(A.13) 
(A.14) 



which proofs equation (1841) . 
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